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1.  INTRODUCTION 


During  the  past  two-and-a-half  years  the  authors  have  sought 

in  this  project  to  develop  an  accurate  and  efficient  tool  for 

predicting  boundary-layer  transition  on  re-entry  vehicle  nose- 

tips.  The  basis  of  our  approach  has  been  a phenomenological 
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turbulence  model  developed  by  Saffman  and  Wilcox.  Unlike 
linear  stability  theory,  the  turbulence  model  approach  includes 
a description  of  the  nonlinear  growth  of  transition  triggering 
disturbances;  the  turbulence  model  approach  thus  provides  a 
description  of  boundary-layer  evolution  from  the  laminar  state, 
through  transition,  and  into  the  turbulence  regime.  In  order 
for  such  a comprehensive  theory  to  remain  tractable,  consider- 
able engineering  judgment  and  approximation  has  been  needed. 
Including  nonlinear  effects,  for  example,  is  achieved  through 
the  introduction  of  several  empirical  closure  coefficients. 

The  notion  that  the  same  values  of  these  coefficients  should 
be  used  for  all  flows  has  been  a central  axiom  in  our  philosophy 
of  developing  the  turbulence-model  transition-prediction  method. 

The  goals  of  our  transition  research  have  actually  been  part  of 
a much  broader  objective.  Specifically,  our  primary  aim  has 
been  to  devise  a set  of  constitutive  equations  suitable  for 
predicting  salient  features  of  practical  engineering  flowfields 
which  are  dominated  by  turbulence  and  transition  phenomena  and 
which  include  complicating  effects  such  as  surface  roughness, 
surface  mass  and  heat  transfer,  boundary-layer  separation, 
streamline  curvature,  etc.  Various  agencies  have  supported 
our  research  efforts  (most  notably  the  NASA  Ames  Research 
Center),  and  as  a result  we  have  made  significant  advances 
toward  achieving  our  basic  objective. 

The  reports  and  papers  listed  in  the  Bibliography  provide 
evidence  of  our  progress  in  treating  the  transition  problem. 
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Using  the  Saf fman-Wilcox  turbulence  model  with  relatively 
straightforward  modifications  to  account  for  low-Reynolds- 
number  effects,  we  have  accurately  predicted  transition  for 
conventional  boundary  layers  and  for  ground  tests  on  re-entry 
vehicle  geometries.  Our  predictions  have  included  most  of 
the  effects  pertinent  to  re-entry  vehicles  including  surface 
roughness,  surface  cooling,  pressure  gradient,  freestream  tur- 
bulence Intensity,  and  freestream  unit  Reynolds  number.  While 
these  results  are  encouraging,  the  number  of  closure  coeffi- 
cients has  increased.  Furthermore,  values  for  some  of  the 
coefficients  have  been  fixed  by  numerical  experimentation, 
i.e.,  by  forcing  agreement  of  theoretical  predictions  with 
experimental  data.  Additionally,  purely  empirical  modifica- 
tions to  the  model  are  required  to  accurately  predict  effects 
of  pressure  gradient  and  surface  heat  transfer  on  boundary- 
layer  transition.  Thus,  even  with  the  model's  recorded  suc- 
cess in  predicting  transition,  there  is  room  for  further 
improvement  of  the  theory. 

An  important  development  in  our  NASA  Ames-sponsored  turbulenc 
research  (Contract  NAS2-8884)  produced  significant  improvement 
in  our  turbulence  model's  accuracy  for  turbulent  boundary 
layers.  Most  significantly,  the  model's  sensitivity  to  ini- 
tial and  boundary  conditions,  a sensitivity  which  has  caused 
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difficulty  in  our  transition  research,  * has  been  reduced 
markedly.  Encouraged  by  this  success,  the  present  study  has 
focused  upon  the  revised  model  equations'  accuracy  for  transi- 
tional flows.  As  will  be  shown  in  the  following  sections,  the 
revised  model  is  superior  to  the  original  model  in  several 
important  respects.  First,  with  no  transition-specific  modifi 
cations,  the  model  accurately  predicts  many  important  features 
of  an  incompressible  flat-plate  boundary  layer  (FPBL)  under- 
going transition,  viz,  velocity  profiles,  skin  friction,  and 
transition  width.  Second,  with  modifications  tied  closely  to 
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linear-stability  theory,  the  model  accurately  predicts  effects 
of  pressure  gradient  and  surface  heat  transfer.  Third,  effects 
of  surface  mass  addition  have  been  incorporated  in  the  theory. 
Perhaps  most  important  however,  while  the  number  of  closure 
coefficients  remains  approximately  the  same  as  in  the  Saffman- 
Wilcox  model,  virtually  all  of  the  coefficients'  values  have 
been  determined  by  methods  more  rigorous  than  numerical 
experimentation . 

Section  2 summarizes  the  model  equations  and  presents  a 
detailed  study  of  viscous  sublayer  structure;  unlike  the  origi- 
nal model,  analysis  of  the  viscous  sublayer  provides  a suit- 
able method  for  treating  low  Reynolds  number  effects.  The 
section  also  discusses  (a)  boundary  conditions  for  surfaces 
with  roughness  and  mass  transfer  and  (b)  modifications  needed 
to  accurately  predict  effects  of  pressure  gradient  and  surface 
mass  addition.  Section  3 presents  boundary-layer  applications 
including  the  incompressible  FPBL,  effects  of  surface  cooling 
on  a supersonic  boundary  layer,  and  nonequilibrium  relaxation 
of  a turbulent  boundary  layer  passing  from  a rough  to  a smooth 
surface.  In  Section  4 , applications  to  blunt  body  flows  are 
presented  including  effects  on  transition  of  surface  rough- 
ness, cooling,  and  mass  addition;  freestream  turbulence 
intensity;  and  freestream  unit  Reynolds  number.  The  concluding 
section  summarizes  results  and  conclusions. 
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2.  THEORETICAL  CONSIDERATIONS 


The  model  equations  are  summarized  in  this  section  including 
established  values  of  all  closure  coefficients.  Subsection  2.1 
presents  physical  meanings  of  turbulence  field  properties  and 
a discussion  of  similarities  and  differences  between  the  new 
model  and  its  predecessor,  the  Saf fman-Wilcox1 ’ ^ model. 
Subsection  2.2  gives  details  of  model-predicted  viscous  sub- 
layer structure,  including  boundary  conditions  for  surfaces 
with  roughness  and  mass  addition.  The  section  concludes  with 
discussion  of  special  modifications  needed  to  obtain  accurate 
predictions  for  effects  on  transition  of  pressure  gradient, 
surface  heat  transfer,  and  surface  roughness. 

2.1  THE  TURBULENCE/TRANSITION  MODEL  EQUATIONS 

Under  the  standard  boundary-layer  approximations,  the  model 
equations  for  two-dimensional  (j=0)  and  axisymmetric  (j=l) 
flows  are 

Mass  Conservation 


|.(pu)  + ^A(rJpv)  0 


(1) 


Momentum  Conservation 


9u  . 9u 

PU37  + pV97  = 


dp  . 3t 
" dx  9y 


(2) 


Energy  Conservation 


3h  . , 9h 

pu37  * »v3y 


d£ 

dx 


9u  3q 

T97  ~ j} 


(3) 
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Turbulent  Mixing  Energy 


PUH  + pv!f  = [a*p ifyl  _ B*w]e  + ^-[(y+a*pe)lf]  (4) 

Turbulent  Dissipation  Rate 

ov!?!‘  {ap|l?1  - [8  + 2c(fy)  ] w|"!  + ^■[(utopc)l7!]  (5> 


where  x and  y are  orthogonal  coordinates  with  x lying  along 

the  body  and  y being  normal  to  the  surface;  r is  the  radial 
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distance  from  the  body  axis.  Mass-averaged  velocity  compo- 
nents in  the  x and  y directions  are  denoted  by  u and  v while 
h is  the  mass  averaged  enthalpy;  p,  p and  p are  mean  density, 
pressure,  and  molecular  viscosity;  t and  q are  the  shear  stress 
and  normal  heat  flux.  The  mass-averaged  turbulent  mixing 
energy,  e,  and  the  mass-averaged  turbulent  dissipation  rate,  w, 
are  needed  to  define  the  eddy  diffusivity,  e,  which  is  given 
by  the  following  equation: 

e = pe/w  (6) 

The  shear  stress  and  heat  flux  are 


T 


(U  + 


„ ! 3u 

pe,5y 


(7) 


q 


_ y_  + _P§_\lk 

\PrL  PrT/3y 


(8) 


where  Pr^  and  PrT  are  laminar  and  turbulent  Prandtl  numbers. 

The  quantity  l is  the  turbulent  length  scale  defined  as 

*.  = peVw  (9) 

The  turbulent  Prandtl  number,  Pr^,  and  the  closure  coefficients 
a, ot*, 8, 8*, o, a#  appearing  in  Equations  (4)  and  (5)  are 
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B = & 
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3_ 

20 


a = 


1 
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B*  = 


r*  = 


100 

_1_ 

2 


Pr  = — 

T 9 


= -i-[l  - (1  - X)  exp  (-ReT/2)] 
= ^[l  - (1  - X)  exp  (-2ReT)] 


(10) 


where  ReT  Is  the  turbulent  Reynolds  number  defined  by 


ReT  = pe^£/y 


(11) 


Specification  of  the  closure  coefficient  X is  deferred  to 
Subsections  2.2  and  2.3. 
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Consistent  with  the  arguments  of  Wilcox  and  Chambers,  the 
turbulent  mixing  energy  is  proportional  to  the  kinetic  energy 
attending  the  fluctuation  of  fluid  particles  normal  to  the 
plane  of  shear.  Letting  v'  denote  the  mass-averaged  fluctuat- 
ing velocity  component  normal  to  the  shear  plane  (under  the 
boundary-layer  approximations,  shear  planes  are  parallel  to 
the  x direction),  the  turbulent  mixing  energy  is  given  by 


9 <pv ’ v ' > 

e = ¥ p 


(12) 


where  p is  the  instantaneous  density. 


The  physical  meaning  of  w has  been  discussed  by  Wilcox  arid 
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Chambers.  For  Incompressible  boundary  layers,  comparison  of 
the  limiting  forms  of  the  model  equations  and  the  exact 
Reynolds  stress  equation  very  close  to  a solid  boundary  shows 
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i r 7* 


that  w/p  is  the  rate  at  which  e is  dissipated  into  heat, 
mean  kinetic  energy,  and  other  fluctuation  modes.  For 
incompressible  flows,  Wilcox  and  Chambers  deduced  that 
w = ( 3p/B* ) < ( 3v ' /3y ) 2 >/<v ' v ' > . In  terms  of  present  notation, 
a suitable  definition  of  w for  compressible  flows  is 


_ 3y  <p(3v'/3y ) 2> 
B * < pv ' v ' > 


(13) 


On  the  one  hand.  Equations  (1-11)  are  very  similar  to  the  model 
equations  developed  by  Saffman  and  Wilcox . 1 ’ ^ The  new 

model  is,  in  fact,  offered  as  the  next-generation  improvement 
over  the  Saf fman-Wilcox  model.  The  most  notable  similarity 
is  the  use  of  the  dissipation-rate,  w,  which  is  analogous  to 
the  frequency  or  pseudovorticity  quantity  first  introduced  by 

Q 

Kolmogorov.  Also , the  high  Reynolds  number  (i.e.,  Re T -*■«•) 
values  of  the  closure  coefficients  defined  in  Equations  (10) 
have  been  determined  by  the  same  arguments  based  on  general 
properties  of  turbulent  flows  used  by  Saffman'*'  and  subse- 

9 

quently  by  Wilcox  and  Alber. 


On  the  other  hand.  Equations  (1-11)  include  several  differences 

from  the  Saffman-Wilcox  model  which  result  in  significant 

improvement  in  accuracy  and  utility.  The  most  important 

difference  is  appearance  of  the  term  proportional  to  (3il/3y)2 

in  Equation  (5).  This  term  has  an  important  impact  on  model- 

predicted  defect-layer  structure  for  an  incompressible  FP3L. 

Physical  meaning  and  origin  of  the  term  have  been  discussed  in 

•3 

detail  by  Wilcox  and  Chambers.  A second  key  difference  is 
In  the  dimensions  of  the  dissipation-rate  quantity,  w.  Its 
dimensions  are  (M/L3t)  which  contrasts  with  the  specific 
pseudovorticity,  ft(L3/Mt),  used  by  Wilcox  and  Alber^  and  the 
pseudovorticity,  w(l/t),  of  Saffman.1  As  shown  by  Wilcox  and 
Traci,  ^ achieving  accuracy  suitable  for  general  engineering 
applications  involving  compressible  flow  is  contingent  upon 
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using  w *v>  poi  rather  than  ui  or  n^w/p.  The  final  noteworthy 


difference  is  in  the  variation  of  a and  a*  with  Re,p.  As 
shown  in  Equations  (10),  a and  a*  approach  their  high  Reynolds 
number  values  for  ReT=0(l).  This  is  physically  more  plausible 
than  a and  a*  achieving  their  Rem  -*■ 00  values  when  Rem=0(10_1) 
as  they  do  in  the  Saffman-Wilcox  model.  * 

In  developing  the  new  model,  Wilcox  and  Chambers  concentrated 
upon  the  high-Reynolds-number  (ReT  + °°)  form  of  Equations  (1-11) 
Values  of  all  the  closure  coefficients  in  the  limit  Re^,  -*■  °° 
have  been  established  by  arguments  based  on  general  properties 
of  turbulent  flows.  The  viscous  modifications,  i.e.,  the  postu 
lated  variation  of  a and  a*  with  Re^  was  devised  in  the  present 
study  by  analyzing  the  viscous  sublayer.  The  next  subsection 
presents  details  of  model-predicted  sublayer  structure  and 
rationale  for  the  viscous  modifications. 

2.2  VISCOUS  SUBLAYER  STRUCTURE 

To  analyze  the  viscous  sublayer,  solutions  are  first  presented 
for  a perfectly  smooth  wall;  the  constant  in  the  law  of  the 
wall,  C,  is  found  to  be  significantly  larger  than  measured 
smooth-wall  values  when  a*  and  a assume  their  high-Reynolds- 
number  values.  Viscous  modifications  for  a and  a*  are  then 
devised  which  yield  a more-reasonable  smooth-wall  value  for  C; 
detailed  comparisons  are  made  between  computed  and  measured 
sublayer  structure.  Next,  a surface  dissipation-rate  boundary 
condition  Is  developed  which  is  suitable  for  turbulent  boundary 
layers  on  rough  surfaces.  Effects  of  mass  injection  on  sub- 
layer structure  are  then  analyzed;  a surface  dissipation-rate 
boundary  condition  is  devised  to  account  for  such  effects  on 
turbulent  boundary  layers.  We  conclude  the  section  by  postu- 
lating a composite  roughness/mass-injection  surface  boundary 
condition  for  the  dissipation  rate. 


8 


2.2.1  Equations  of  Motion 

In  the  viscous  sublayer  of  an  incompressible  turbulent  boundary 
layer,  convective  terms  are  negligible  so  that  the  equations 
of  motion  become 

(v  + e)^  = u*  (14) 


[6  + 20  (a!)]  5) 


w2  + -r-  (v+ae) 


1- 


where  v is  kinematic  viscosity  and  u = A /p  is  friction  veloc- 

t w 

Ity.  Five  boundary  conditions  must  be  specified  for  this  fifth- 
order  set  of  ordinary  differential  equations.  Two  boundary 
conditions  follow  from  asymptotic  behavior  of  solutions  to  the 
model  equations  as  u^y/v-*-00.  Noting  that  turbulent-energy 
diffusion  is  negligible  in  this  limit,  there  follows 


u2/a* 

X 00 


pu  /a#icy  as  u y/v  -*•  00 

OO  v 


where  a*  = 3/10  is  the  limiting  value  of  a*  for  Re^-*-00.  Note 
that  Equations  (17)  are  consistent  with  the  law  of  the  wall, 
i . e . , 

— = I log  lid  + c (: 

U K 6 V 

T 

where  k = 0.41  is  Karman's  constant. 


Two  surface  boundary  conditions  follow  from  the  no-slip  bound- 
ary condition  which  implies  that  u and  e vanish  at  y=0. 

1 2 

Following  Saffman  and  Saffman  and  Wilcox,  we  impose  a sur- 
face boundary  condition  on  the  dissipation  rate  so  that 


u = e = 0 


H ■„(¥) 


at  u^y/v  = 0 


where  SR  is  a universal  function  of  the  roughness  height  of 
the  wall,  k.  As  shown  by  Saffman  and  Wilcox,  the  constant  in 
the  law  of  the  wall,  C,  depends  upon  SR.  Since  experimental 
observations  show  that  C also  depends  upon  u^k/v,  comparison 
of  computed  and  measured  values  of  C establishes  the  dependence 
of  SD  upon  u k/v. 

n T 


2.2.2  Smooth-Wall  Sublayer  Structure 

To  examine  model-predicted  sublayer  structure,  we  first  con- 
sider the  case  of  a perfectly  smooth  wall,  i.e.,  SR  -*■«>.  One 
of  the  most  significant  results  obtained  by  Wilcox  and  Saffman 
is  that,  with  no  viscous  modifications  to  their  model  other 
than  including  molecular  diffusion,  the  constant  C -*  5 • 5 for  a 
perfectly  smooth  wall.  By  comparison,  measurements  indicate 
that  C should  lie  somewhere  between  5 and  6.  Note  that 
Saffman  and  Wilcox  deduced  C->5-7  for  perfectly  smooth  walls. 
However,  in  this  study,  Saf fman-Wilcox-model  computations  were 
done  using  more  accurate  methods  than  those  employed  by  Saffman 
and  Wilcox.  Our  computations  Indicate  C-+5-5  in  the  limit  of 
a perfectly  smooth  wall. 


For  a given  value  of  SD  the  corresponding  value  of  C is  deter- 

n 

mined  by  first  solving  the  two-point  boundary-value  problem 
defined  by  Equations  (1^4-17,  19),  and  by  then  evaluating  the 
following  limit: 


C 


lim 
y+  -*•  °° 


(20) 


where  y+  = uxy/v . 


The  boundary  value  problem  was  solved  in  the  present  study 
with  an  implicit,  second-order-accurate,  time-marching,  finite 
difference  method.  Initially,  computations  were  performed 
with  a*  = a*  and  a = a^.  Extrapolation  of  the  numerical  results 


10 


for  large-but-f inite  Sp  indicated  that,  with  the  new  model, 
the  limiting  value  of  C for  a perfectly  smooth  wall  is  approx- 
imately 7-0.  Such  a large  limiting  value  of  C is  unacceptable 
for  general  engineering  applications,  thus  prompting  develop- 
ment of  further  viscous  modifications  to  the  model  equations. 


2.2.3  Viscous  Modifications 


Suitable  viscous  modifications  are  developed  by  noting  that, 

5 

as  argued  by  Wilcox  and  Chambers,  when  Re^ ^ 1 net  production 
of  turbulent  energy  is  reduced  relative  to  the  Rem  >>  1 situa- 

5 1 

tion.  This,  leads  Wilcox  and  Chambers,  in  a prior  study  of 
boundary-layer  transition,  to  propose  a viscous  modification 
to  the  Saffman-Wilcox  model  which  reduces  the  production  of 
turbulent  mixing  energy  for  small  Re^.  As  a generalization  of 
Wilcox's  and  Chambers'  viscous  modification,  we  propose  that 
a and  a*  depend  upon  ReT  as  follows : 


a*  = a * [l  - (1  - X)  exp  (-ReT/Re)] 
a = [l  - (1  - X)  exp  (-ReT/Rw)] 


(21) 


Equations  (21)  contain  three  closure  coefficients,  i.e.,  X,  R , 
and  R . The  value  of  X can  be  determined  by  demanding  that 
the  model  equations  predict  that  in  a Blasius  boundary  layer 
turbulent  fluctuations  are  damped  for  Reynolds  numbers  below 
the  linear-stability-theory  minimum-critical  Reynolds  number, 
Rec.  Having  mixing-energy  production,  a*|9u/3y|e,  less  than 
mixing  energy  dissipation,  g*we/p,  insures  such  damping. 

Using  the  Blasius  velocity  profile  and  the  smooth  wall  w pro- 
file  (i.e.,  w = 20p/3y  — see  Wilcox  and  Wilcox  and  Chambers  ) , 

the  maximum  plate-length  Reynolds  number,  Re-',  at  which  dis- 
sipation is  greater  than  or  equal  to  production  throughout 
the  boundary  layer  is 

Re-  = 750/X2  (22) 


M- 
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the  accepted 


The  Reynolds  number  Re^  will  be  equal  to  9*10“, 
value  of  Rec,  provided 

X = 1/11  (23) 

Note  that  in  order  to  accurately  predict  effects  of  pressure 
gradient  and  surface  heat  transfer.  Equation  (23)  must  be 
modified.  Further  details  are  given  in  Subsection  2.3. 

Values  of  Rg  and  were  chosen  to  achieve  optimum  agreement 
between  computed  and  measured  sublayer  structure.  As  the 
first  step  in  determining  the  optimum  (R  ,Rw)  pair,  computa- 
tions were  performed  for  various  values  of  Rg  and  Rw  to  deter- 
mine pairs  for  which  the  smooth-wall  value  of  C is  within  the 
measured  range.  As  shown  in  Figure  1,  a locus  of  values  (RgjR^ 
exists  which  yields  C=5-5  for  smooth  walls. 

The  next  step  was  to  examine  predicted  sublayer  profiles  and  to 
compare  with  experimental  sublayer  data  of  Laufer.^  The  com- 
parisons indicated  that  closest  agreement  with  Laufer's  data 
is  obtained  with 


R. 


R 


w 


1/2 

2 


(24) 


Figure  2 compares  computed  and  measured  velocity  and  shear- 
stress  profiles  and  illustrates  the  close  agreement  between 
the  computed  and  measured  turbulent-energy  production/ 
dissipation  balance. 

As  noted  in  Subsection  2.1,  a and  a*  approach  their  fully 
turbulent  (Re,p>>  1)  values  for  Re,p  = 0(1)  while,  in  the  Saffman- 
Wilcox  model,  a = an  and  a*  = a*  for  Re,p=10_1.  Intuitively, 
it  seems  unrealistic  that  a flow  should  begin  to  exhibit  char- 
acteristics of  fully  developed  turbulence  before  the  eddy 
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diffusivity  even  becomes  comparable  to  the  molecular  diffu- 

sivity.  The  viscous  modifications  developed  here  thus  appear 

4 

more  realistic  than  those  devised  by  Wilcox  which,  in  effect, 

have  R = R = 1/10.  Wilcox  was  forced  to  use  very  small  values 
e w 

of  R and  R because  larger  values  yielded  too  small  of  a 
0 w 

smooth-wall  C.  This  is  unsurprising  in  light  of  Figure  1. 
Conceivably,  a similar  locus  of  (Re>Rw)  pairs  exists  for  the 
Saf fman-Wilcox  model,  although  no  attempts  have  been  made  to 
find  such  pairs.  The  fact  that  the  locus  passes  through  the 
origin  for  the  Saf fman-Wilcox  model  is  somewhat  fortuitous. 

2.2.4  Rough-Surface  Boundary  Condition 

Having  achieved  acceptable  smooth-wall  sublayer  structure,  we 

now  proceed  to  determine  the  dependence  of  SR  upon  u^k/v.  As 

noted  earlier  in  this  section,  C depends  upon  SR  which,  in 

turn,  is  a function  of  surface  roughness.  Figure  3 shows  the 

predicted  variation  of  C with  SD;  results  are  similar  to  those 

n 12 

obtained  by  Saffman  and  Wilcox.  Sufficient  experimental  data 

are  available  for  effects  of  "sand-grain"  roughness  on  turbu- 
lent boundary-layer  velocity  profiles  to  define  C uniquely  as 
a function  of 


k+  = uTk/v  (25) 

Comparing  the  measured  variation  of  C with  k+  to  the  computed 
variation  of  C with  SR  leads  to  a correlation  for  SR  as  a 
function  of  k+;  Figure  3 shows  the  correlation.  An  accurate 
analytical  fit  is  given  by 

SR  = (36/k+)2  + (8/k+)Ss  (26) 

2.2.5  Effects  of  Mass  Injection^ 

For  boundary  layers  with  surface  mass  injection  (blowing), 
the  introduction  of  an  additional  velocity  scale  (v  = normal 

^ Results  discussed  in  this  sub-subsection  were  obtained  under 
sponsorship  of  the  NASA  Langley  Research  Center  (Contract  NAS1- 
13974). 
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Figure  3-  Basis  of  correlation  between  surface  dissipation-rate  boundary 
condition  and  surface  roughness. 


flow  velocity  at  the  surface)  suggests  that  some  modification 
of  Equation  (19)  may  be  required  for  flows  with  blowing. 
Andersen,  et  al,1^  provide  further  evidence  that  the  dissipa- 
tion-rate boundary  condition  must  be  revised  when  blowing  is 
present  by  showing,  from  correlation  of  their  experimental 
data,  that  the  law  of  the  wall  assumes  the  following  modified 
form: 


u 


+ C 


(27) 


The  modified  "constants"  < and  C are  related  to  < , C,  and  the 
blowing  velocity  as  follows: 


< = K/(l  + 7.7  vw/uT) 

C = C - 50(vw/ut)  + 72(vw/ut)2 


(28) 


Since  the  computations  above  show  that  C is  strongly  affected 
by  the  value  of  S^,  Andersen's  data  indicate  that  a modification 
to  the  model  equation  boundary  condition  is  required  for  com- 
puting blown  boundary  layers.  Since  mass  injection  is  often 
pertinent  for  re-entry  vehicles,  we  thus  consider  effects  of 
blowing  on  model-predicted  sublayer  structure. 


When  blowing  is  present  the  sublayer  equations  are 

(v+£)57  = ux  + vwu  (29) 

- H£l  - MS}.  ♦&[(*♦...>£]  (30, 

% r - {°i§i  - [»«•($]*}■•  * tl  <3i> 

To  establish  boundary  conditions  for  this  fifth-order  system, 
we  assume  the  effluent  gas  is  free  of  turbulent  fluctuations 
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so  that  e vanishes  at  the  surface.  The  horizontal  velocity,  u, 
vanishes  by  virtue  of  the  no-slip  condition,  and  in  analogy  to 
Equation  (19)  we  write 


u = e = 0 


w 


uTy/v 


0 


(32) 


where  is  a universal  function  of  v, ,/u  . Examination  of 
B w x 

asymptotic  solution  behavior  yields  two  more  boundary  conditions 
valid  as  u y/v -*•“>.  As  in  the  no-blowing  case,  the  turbulent 
diffusion  term  in  the  mixing-energy  equation  is  negligible  for 
uTy/v  -*■  °°,  and  there  follows 


u* 


v 
u u 


1 + 


T J 


w 


p du 
a*  dy 


as 


uTy/v 


(33) 


In  addition  to  establishing  boundary  conditions,  further  exami- 
nation of  the  wall-layer  solution  (i.e.,  the  limiting  solution 
for  u.j.y/v-*00)  demonstrates  direct  correspondence  between  model- 
predicted  and  measured  effects  of  blowing  on  a turbulent  bound- 
ary layer.  Specifically,  in  the  limit  of  weak  blowing  (v  /u  <<  1 
expanding  in  powers  of  v /u  shows  that  the  velocity  obeys  a modi 
fled  law  of  the  wall  similar  to  Equation  (27).  The  effective 
Karman  constant  is  predicted  to  be 


ic  = k/(1+5vw/ut  ) 

where 


4kC-1  . 
" 8< 


uty 


(34) 


(35) 


The  functional  dependence  in  Equation  (34)  is  similar  to  that 
quoted  by  Andersen  [Equation  (28)].  Table  1 shows  the  varia- 
tion of  E with  u^y/v  for  < = 0.4l  and  C®5.5*  As  shown,  the 
predicted  value  of  E is  reasonably  close  to  Andersen's  value 
of  7-7,  particularly  at  the  larger  values  of  uTy/v. 


Table  1.  Variation  of  Z with  u^y/v. 


uTy/v 

■z 

10 

3.84 

100 

5.^5 

500 

6.23 

1000 

6.65 

Again  using  the  implicit,  time-marching  numerical  method,  sub- 
layer calculations  were  performed  for  several  blowing  rates 
ranging  from  vw/uT  = 0 to  vw/uT  = 0.7393.  The  value  of  Sg  was 

varied  for  each  blowing  rate  until  close  agreement  with  the 

13 

velocity  profile  data  of  Andersen,  et  al , was  obtained 

(Figure  4 ) . As  expected  from  the  wall-layer  analysis,  slopes 

of  the  various  profiles  (l.e.,  ic_1)  are  accurately  predicted. 

Note  that  this  means  that,  similar  to  the  case  of  surface 

roughness,  the  dissipation-rate  boundary  condition  primarily 

determines  the  variation  of  C with  v/u  . Figure  4 also  pre- 

w x 

sents  a correlation  of  Sg  with  vw/uT ; an  accurate  analytical 
representation  of  the  correlation  is 


S 


B 


6 


1 + <VUT> 


v /u 

W T 


> 0 


(36) 


Although  effects  of  suction  (v/u  <0)  have  not  been  considered 

w x j4 

here,  computations  performed  by  Wilcox  imply  that  w is 
unaffected  by  suction.  Therefore,  Equation  (36)  should  only 
be  used  for  v/u_  > 0. 

W T 


2.2.6  Composite  Roughness/Blowing  Boundary  Condition 

We  conclude  this  section  on  model-predicted  sublayer  structure 
by  postulating  a dissipation-rate  boundary  condition  suitable 
for  surfaces  with  both  roughness  and  blowing.  Physical  con- 
siderations must  guide  development  of  such  a boundary  condition. 


C*  ^ 
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COMPARISON  OF  COMPUTED  AND  MEASURED  CORRELATION  OF  SR  WITH  BLOWING 
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Figure  Basis  of  correlation  between  surface  dissipation-rate  boundary 
condition  and  surface  mass  addition  rate. 


For  example,  on  extremely  rough  surfaces  (k+-*-°°)  we  expect 
the  kinetic  energy  of  the  injected  fluid  to  be  small  compared 
to  the  energy  attending  local  flow  separations  and  turbulent 
fluctuations  generated  by  the  roughness  elements.  Effects  of 
finite  blowing  rates  are  therefore  unimportant  in  this  limit. 
Similarly,  for  large  blowing  rates  (v  /u  -*•<*>),  finite  rough- 

W T 

ness  heights  have  a negligible  effect.  Additionally,  the 
composite  boundary  condition  must  reduce  to  (a)  Equations  (19) 
and  (26)  for  v /u  = 0 and  (b)  Equations  (32)  and  (36)  for 
k+=0.  The  following  boundary  condition  satisfies  all  of 
these  constraints: 

pux 

w = ^S  at  y = 0 (37) 

oo 

where 

s = (^  + 4)’1  <38) 

with  SR  and  Sfi  given  by  Equations  (26)  and  (36). 

At  this  point,  the  model  has  been  developed  with  virtually  no 
transition-specific  considerations.  That  is,  with  the  excep- 
tion of  the  argument  establishing  the  value  of  X [Equation  (23)], 
values  of  all  closure  coefficients  have  been  fixed  by  arguments 
based  on  properties  of  fully  turbulent  flows.  As  will  be  shown 
in  Section  3»  the  model,  with  no  further  modifications,  pro- 
vides an  accurate  description  of  many  subtle  aspects  of  incom- 
pressible, zero-pressure-gradient  boundary-layer  transition 
(Section  3).  However,  similar  to  the  Saffman-Wilcox  formula- 
tion, the  new  model  requires  modifications  to  the  closure 
coefficient  X in  order  to  accurately  predict  effects  on  tran- 
sition of  pressure  gradient  and  surface  heat  transfer;  addi- 
tionally, the  rough-surface  boundary  conditions  must  be  revised 
for  transition  applications.  Suitable  modifications  to  X, 
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closely  tied  to  linear  stability  theory,  have  been  devised 
under  joint  sponsorship  (Contract  N00024-76-C-7070 ) of  the 
Defense  Advanced  Research  Projects  Agency  (DARPA)  and  the 
Naval  Sea  Systems  Command  (NAVSEA).  Because  of  the  pertinence 
of  the  modifications  to  the  present  study,  the  following  sub- 
section presents  complete  details  before  proceeding  to 
applications  in  Sections  3 and 

2.3  TRANSITION  MODIFICATIONS 

In  the  first  part  of  this  subsection  we  focus  on  the  transi- 
tion modifications  to  the  closure  coefficient  A.  The  second 
part  presents  rough-surface  boundary-condition  revision  needed 
for  transition  computations. 

2.3-1  Pressure  Gradient  and  Heat  Transfer 

As  noted  in  the  preceding  section,  the  value  of  A has  been 
fixed  by  demanding  that  the  linear-stability  minimum-critical 
Reynolds  number,  Rec,  for  the  Blasius  boundary  layer  match 
the  corresponding  model-equation  neutral-stability  Reynolds 
number.  Re-'.  Demanding  that  Re  = Re-'  for  the  Blasius  boundary 
layer  yields  the  value  of  A given  in  Equation  (23).  The  model 
equations  reasonably  can  be  expected  to  apply  to  transitional 
flows  which  are  insensitive  to  spectral  effects.  That  is,  the 
various  constants  in  the  model  equations  are  essentially  cor- 
relation coefficients  which  have  been  integrated  over  the 
turbulent  spectrum.  Hence,  if  the  stability  diagram  shows 
that  a wide  range  of  wave  numbers,  a,  undergo  amplification, 
the  spectrum  will  more  closely  resemble  a fully-turbulent 
spectrum  than  if  only  a small  range  of  wave  numbers  are 
unstable.  For  example,  the  stability  diagram  for  a boundary 
layer  subjected  to  a pressure  gradient  is  shown  in  Figure  5. 
For  adverse  pressure  gradient,  a finite  range  of  wave  numbers 
are  unstable  at  all  Reynolds  numbers  in  excess  of  ReXc  (note 
that  6*  is  displacement  thickness).  On  the  basis  of  the  dis- 
cussion above,  the  model  would  be  expected  to  accurately 
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Figure  5.  Neutral  stability  curves  for 

laminar  boundary  layers  with  pres 
sure  gradient  and  suction. 


predict  the  destabilizing  effect  of  adverse  pressure  gradient. 
In  contrast,  the  stability  diagram  becomes  thinner  with 
increasing  favorable  pressure  gradient  so  that  spectral 
effects  become  increasingly  important,  particularly  for  small 
T'  which  yields  transition  at  large  values  of  Re^#;  the  model 
hence  would  be  expected  to  fare  poorly  for  transitional  bound- 
ary layers  with  favorable  gradients  (and  low  freestream 
disturbances ) . 

c; 

As  shown  by  Wilcox  and  Chambers,  the  original  version  of  the 
model  behaves  just  as  the  above  discussion  indicates.  Excel- 
lent agreement  between  theory  and  experiment  is  obtained  for 
adverse  gradients  while,  for  low  freestream  turbulence  inten- 
sities, the  model  fails  to  predict  the  strong  stabilizing 

effect  of  favorable  gradient.  To  remove  this  deficiency, 

5 

Wilcox  and  Chambers^  introduced  an  empirical  modification  to 
Equation  (23).  While  reasonably  good  agreement  with  measure- 
ments resulted,  the  modification  lacked  rigor. 


In  the  NAVSEA  study,  a better  approach  to  modifying  X has  been 
found.  That  is,  the  requirement 

Rec  = Re-  (39) 


has  been  extended  to  include  favorable  pressure  gradients. 
Hence,  for  small  freestream  turbulence  intensity  we  expect  to 
have 

X = yi  f ( A)  as  T'  -*■  0 (40) 


where  A is  the  modified  Pohlhausen  pressure  gradient  parameter 

P « ft*  P 


e = _ ^e  0^/32u\ 

P, . V ~33T  ” p U \3y i ) 

K w e we  w 


(41) 


The  function  f(A)  must  be  determined  by  equating  Re„  and  Re-. 

v A 

Figure  6 presents  results  based  on  the  Pohlhausen  profiles!**; 
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Figure  6.  The  stability  function  f(A). 
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a good  fit  to  the  data  indicate  the  variation  of  A with  A is 
hence 

f ( A)  = ^ + |jexp  [-40AH(A)]  (42) 

where  H(A)  is  the  Heaviside  stepfunction . The  limiting  value 
of  f(A)  as  A + 0 has  been  obtained  by  a similar  analysis  of  the 
asymptotic  laminar  profile  for  a uniformly -sucked  FPBL. 


I 

i 


Turning  to  heat  transfer,  an  additional  modification  to  A is 
needed.  Following  Wilcox  and  Chambers, ^ for  incompressible 
aerodynamic  boundary  layers  with  heat  transfer,  the  neutral 
stability  Reynolds  number  based  on  wall  conditions  is  given 


by 


750 

"A3" 


(43) 
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According  to  linear  stability  theory,  the  corresponding 

minimum-critical  Reynolds  number  ReP  varies  as  follows: 

cw 

Rec  ^ (W"7  (41° 

w 

Hence,  since  y ^ T0*7  for  air,  we  have 

Recw  * (^w/yer5  (45) 

We  thus  postulate  that  for  flows  with  heat  transfer,  A be 
written  as  follows: 


(46) 


Finally,  note  that  for  high  freestream  turbulence  intensities, 
transition  is  unlikely  to  be  sensitive  to  spectral  effects 
regardless  of  the  stability  diagram.  That  is,  typical  high 
intensity  freestream  turbulence  will  have  fluctuations  at 


all  frequencies  (wave  numbers).  Therefore,  the  modifications 

proposed  in  Equation  (46)  is  strictly  valid  only  at  T'+O. 

Hence,  to  complete  the  formulation , we  introduce  an  exponential 

dependence  upon  T'  , the  maximum  disturbance  in  the  boundary 

max 

layer,  so  that  the  completed  formulation  is  as  follows: 


TRANSITION  MODIFICATION 


- 1 + tf  CA)-i:  exp 
r(A>  = — + ||exp  [-“OAHU)] 


where 
A = 


T’ 

max 


100 


max 


/U 


e 


(47) 


Figures  7 and  8 illustrate  how  accurately  the  model  predicts 

effects  on  transition  of  pressure  gradient  and  heat  transfer. 

Figure  7 compares  predicted  and  measured  effects  of  pressure 

gradient  on  FPBL  transition.  As  shown  in  the  figure,  the 

Crabtree1^  data  are  closely  simulated  by  the  computations 

with  T'  = .01%  and  .03%.  Since  the  Crabtree  data  are  for 

flight  tests  (closed  circles)  and  quiet  wind  tunnels  (open 

circles),  the  value  of  T*  for  the  experiments  would  be 

expected  to  fall  in  the  range  of  .01-. 03%.  The  agreement 

between  theory  and  experiment  hence  is  very  good  for  low 

17 

intensities.  Additionally,  the  Feindt  high-intensity  data 

(T1  =1.25%)  are  closely  matched  by  the  computed  curve  with 

T'  =1.25%.  Figure  8 compares  computed  and  measured  transition 

Reynolds  number  for  low-speed  aerodynamic  boundary  layers. 

1 8 

The  data  of  Zysina-Mc! ozhen  and  Kuznetsova  were  taken  in 
relatively  noisy  environments  so  that  they  correspond  to  high 
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Figure  7.  Comparison  of  computed  and  measured 

effects  of  pressure  gradient  on  boundary 
layer  transition. 


intensity  freestream  conditions.  The  fact  that  the  curve 
computed  with  T'  = 1.25?  Is  close  to  the  data  is  hence  very 
encouraging.  For  the  low  Intensity  computations  (T*  = 0.03/0, 
the  predicted  stabilizing  effect  of  cooling  is  much  stronger 
than  for  the  T'  =1.25?  computations.  As  expected,  the  varia- 
tion of^ReXt^  is  approaching  the  linear  stability  predictions 
of  Mack. 


2.3.2  Rough-Surface  Boundary  Conditions 

As  with  the  Saffman-Wilcox  model,  the  rough-surface  boundary 

conditions  for  e and  w must  be  revised  for  transitional  flows. 

The  primary  reason  for  the  modifications  is  to  facilitate  use 

of  a boundary  layer  computation  for  predicting  transition;  use 

of  Equations  (19)  would  require  an  elliptic  integration  method 

to  account  for  local  flow  separation  between  roughness  elements. 

Complete  discussion  of  this  point  is  given  by  Wilcox  and 
5 

Chamb  ers . 

4 

Similar  to  the  approach  taken  by  Wilcox  and  by  Wilcox  and 
5 

Chambers,  we  postulate  that  Equation  (19)  be  replaced  by 


where  SD  is  given  by  Equation  (26).  The  functions  N and  N 
n e w 

have  been  found  by  numerical  experimentation;  their  postulated 
dependence  upon  cf./(k/0)  is 


k+  < 5 
k+  > 5 


(49) 
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f 


k/0  < 275  c , 


N = 
w 


r-275  cf-|6 

L k/0  J 


k/0  > 275  c . 


(50) 


Equations  (48-50)  closely  resemble  the  boundary  condition 
revisions  presented  by  Wilcox  and  Chambers.^  The  most  signif- 
icant difference  is  in  Equation  (49),  where  e is  forced  to 
vanish  for  k+<  5-  Having  such  a cutoff  is  equivalent  to  stat- 
ing that  below  a critical  Reynolds  number,  vortex  shedding 
behind  roughness  elements  ceases.  The  cutoff  value  of  5 has 

been  established  by  the  following  argument.  First,  experi- 

19 

mental  measurements  indicate  that  for  flow  past  a cylinder, 
vortex  shedding  occurs  only  when  Reynolds  number  based  on 
cylinder  radius  exceeds  about  25.  Second,  the  corresponding 
Reynolds  number  for  a roughness  element  is  u^k/v  where  u^  is 
a characteristic  velocity.  Third,  assuming  the  velocity 
varies  linearly  close  to  the  wall,  then  u.  ^u2k/v;  conse- 
quently,  u^k/v = (k+) 2 where  k+  is  defined  in  Equation  (25). 
Finally,  requiring  a cessation  of  vortex  shedding  for 
u^k/v  < 25  implies  no  vortex  shedding  for  k+<  5- 


Figure  9 shows  computed  effects  of  surface  roughness  or.  FPBL 

17 

transition;  experimental  data  of  Feindt  are  shown  for  com- 
parison. Freestream  intensity  has  only  a slight  effect  on 
transition  Reynolds  number  based  on  plate  length,  ReXfc , for 
roughness-height  Reynolds  numbers  in  excess  of  300.  Thus, 
consistent  with  qualitative  observations,  the  revised  model 
predicts  existence  of  a roughness  dominated  regime,  although 
substantiating  data  are  unavailable  for  determining  the 
minimum  value  of  Re^  at  which  transition  becomes  roughness 
dominated.  Also,  again  consistent  with  measurements  and 
qualitative  observations,  roughness  has  virtually  no  effect 
on  FPBL  transition  for  Re^$120. 
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Figure  9 • Effect  of  surface  roughness  on  flat-plate 
boundary-layer  transition. 


The  turbulence/transition  model  is  now  completely  formulated. 
In  summary  (a)  -the  completed  set  of  model  equations  is  given 
by  Equations  (1-11),  (b)  the  closure  coefficient  X is  defined 
in  Equation  (U7),  and  (c)  surface  boundary  conditions  for  e 
and  w are 


where 


(51) 


(52) 


with  SR  and  SR  given  in  Equations  (26)  and  (36),  respectively; 
N and  N are  defined  by  Equations  (^9)  and  (50).  All  of  the 

6 W 

applications  in  the  following  two  sections  have  been  done  with 
no  adjustment  of  the  various  closure  coefficients. 


3.  BOUNDARY  LAYER  APPLICATIONS 


To  test  the  new  model,  we  first  consider  carefully-documented 
experiments  for  conventional  boundary-layers.  This  section 
presents  results  of  three  applications.  First,  we  apply  the 
model  to  an  incompressible  FPBL  including  comparisons  of 
computed  and  measured  transitional  velocity  profiles,  transi- 
tion width,  and  transition  sensitivity  to  freestream  turbulence 
intensity  and  scale.  Then,  we  simulate  the  stabilizing  effect 
of  surface  cooling  on  a supersonic  boundary  layer.  Finally, 
we  predict  the  nonequilibrium  adjustment  of  a rough-wall 
turbulent  boundary  layer  which  is  suddenly  exposed  to  a 
smooth  surface. 

3.1  THE  INCOMPRESSIBLE  FLAT  PLATE  BOUNDARY  LAYER 

As  the  first  step  in  testing  the  new  model,  we  analyze  various 
aspects  of  model-predicted  incompressible  FPBL  transition.  As 
shown  in  Figure  10,  starting  from  laminar  flow  at  the  plate 
leading  edge,  the  model  predicts  that  the  skin  friction  ini- 
tially matches  the  Blasius  value.  Then,  depending  upon  the 
freestream  turbulence  intensity,  T' , defined  by 

T'  = 100  ^ !%/Ue  (53) 

the  skin  friction  rapidly  increases  at  a critical  Reynolds 
number,  ReXt,  and  asymptotically  approaches  the  equilibrium 
turbulent  value.  Predicted  variation  of  c^.  closely  resembles 
that  observed  when  a boundary  layer  undergoes  transition  to 
turbulence . 

For  example.  Figure  11  shows  that,  consistent  with  measure- 
20 

ments,  momentum- thickness  Reynolds  number  at  transition, 

Reg  , varies  almost  linearly  with  T'  for  low-intensity  free- 

V/ 

stream  turbulence  (i.e.,  T'  less  than  1%) . (The  criterion 
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Figure  10.  Computed  variation  of  skin  friction  with  Rex  for  transitional  boundary 
layers;  transition  region  lies  between  points  at  which  minimum  and 
maximum  c-,  occur. 


Figure  11.  Effect  of  freestream  turbulence  intensity  on  transition 
location  for  a flat-plate  boundary  layer. 


used  to  define  transition  is  the  point  at  which  c^.  achieves 
a minimum. ) Note  that  turbulence  scale  has  an  effect  on 
ReQ  , particularly  for  T'  >1?.  Computations  have  been  per- 
formed with  various  ratios  of  i.e//a£  to  6 ; having 
given  by 

= .004  ; T'  < 1 (54) 

/a  *'  6 

CO 


most  nearly  matches  the  low-intensity  values  of  Reg^_.  Some- 
what larger  values  of  SLe/  6 are  needed  to  match  the  high- 
intensity  data;  for  values  of  T'  in  excess  of  1% , excellent 
agreement  between  computed  and  measured  Reg^  obtained 
with  &e//o J <S  given  by 


.004  T*  ; T"  > 1 


(55) 


Intuitively,  we  expect  that  should  increase  with  T'  for 
high-lnensity  turbulence  since,  in  the  case  of  fully  turbulent 
boundary  layers,  values  of  l / /a**  6 generally  are  an  order-of- 

C 00 

magnitude  greater  than  that  given  by  Equation  (5^).  The  vari- 
ation of  le  with  T*  given  by  Equation  (55)  is  quite  plausible 
since  the  peak  local  intensity  in  a fully  turbulent  incompres- 
sible FPBL  is  of  the  order  of  10  to  12%  which,  from  Equa- 
tion (55)  indicates  is  of  the  order  .04  to  .05-  For 

G 00  

turbulent  boundary  layers,  i / /a£ 6 is  typically  .09. 

G w 


In  Figure  12,  predicted  width  of  the  transition  region  is 

21 

compared  with  measured  width.  Transition  width,  Axt , is 
defined  as  the  distance  between  minimum  and  maximum  skin- 
friction  p'.ints.  Using  this  definition  for  Axfc,  Reynolds 
number  based  on  Ax.  was  computed  for  Rev  ranging  from 

Z 

5.0*101'  to  4.4*106.  As  shown,  the  computed  curve  falls 
within  experimental  data  scatter. 
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Figure  13  presents  a direct  comparison  between  computed  and 
22 

measured  transitional  total  skin  friction,  Cp,  defined  by 


C 


F 


D 


(56) 


where  D is  drag  per  unit  area  and  x is  distance  from  the  plate' 

leading  edge.  Consistent  with  the  data,  the  boundary  layer  is 

predicted  to  have  a long  memory  of  upstream  flow  conditions. 

That  is,  while  transition  occurs  at  Rex=4*10s,  the  predicted 

Cp  doesn't  begin  its  asymptote  to  the  fully  turbulent  value 

until  4*10®,  a full  decade  higher.  The  long-persisting  memory 

of  transition  is  also  indicated  in  Figure  10  where  the  local 

skin-friction  coefficient  first  overshoots  and  then  gradually 

approaches  the  fully  turbulent  value.  Such  an  overshoot  in 

21 

Cp  is  observed,  particularly  in  supersonic  boundary  layers  ; 
the  Cp  overshoot  is  believed  to  be  associated  with  the  mea- 
sured higher-than-expected  recovery  factors  near  transition. 

21 

Figure  14  compares  the  predicted  and  measured  magnitude  of 
the  overshoot  defined  by 


(57) 


22 

where  cp  is  given  by  the  Karman-Schoenherr  correlation  for 
turbulent  FPBL's.  As  can  be  seen  from  the  figure,  model  pre- 
dictions are  close  to  the  experimental  data.  A good  fit  to 
the  numerical  results  is 


Ac„ 

— ~ = .063  log  Re  - B 

Cf  xt 

o 


(58) 


where  B=  0.734  when  the  transition  point  is  defined  as  the 
point  where  Cp  achieves  a minimum  and  B=  0.720  when  the  tran- 
sition point  is  defined  as  the  point  where  Cp  first  differs 
from  the  laminar  value  by  more  than  0.55&. 
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Figure  13.  Comparison  of  computed  and  measured  total  skin  friction  for 
transitional  flat-plate  boundary  layer. 
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Figure  14.  Comparison  of  computed  and  measured 
skin  friction  overshoot  through 
transition  for  a flat-plate  boundary 
lay 
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Figure  15  exhibits  computed  and  measured  velocity  profiles 
through  transition  for  T'  = .03%.  Comparison  of  velocity  pro- 
files at  x=5-75  ft  indicates  the  computed  boundary  layer 
goes  turbulent  a bit  faster  than  measured.  Proceeding  down- 
stream, however,  computed  and  measured  velocity  profiles 
show  decreasing  differences  until,  beyond  x=6.75  ft,  differ- 
ences are  less  than  5%. 

An  interesting  point  about  these  results  is  that  because 

(32u/3y2)  vanishes  for  incompressible  FPBL  flow,  the  transi- 

w 

tion  modification  defined  in  Equation  (47)  has  no  effect. 

Hence,  the  turbulence  model  as  stated  in  Subsection  2.1  applies, 
with  no  alteration,  to  the  incompressib_e  FPBL. 

3.2  COMPRESSIBLE  BOUNDARY  LAYER  WITH  HEAT  TRANSFER 

Our  second  application  is  for  the  effects  of  surface  cooling 

on  a smooth-wall,  Mach  2.7  boundary  layer.  Figure  16  compares 

computed  effects  of  cooling  with  the  experimental  data  of 

25 

van  Driest  and  Boison  , T and  T denote  wall  temperature 

w aw  ^ 

and  adiabatic  wall  temperature,  respectively.  Computations 
have  been  done  for  low  (T'  = 0.01%)  and  high  (T'  = 1.25%)  free- 
stream  turbulence  intensities.  Consistent  with  measurements, 
cooling  has  a strong  stabilizing  effect,  particularly  for  low 
freestream  turbulence  level.  Note  that,  consistent  with  the 
van  Driest-Boison  data,  the  effectiveness  of  cooling  is  greatly 
reduced  when  the  freestream  turbulence  is  very  intense. 

3-3  BOUNDARY  LAYER  SUBJECT  TO  A STEP  CHANGE  IN 
ROUGHNESS 

To  test  the  model's  ability  to  predict  effects  of  surface 
roughness,  we  apply  the  model  to  the  development  of  a turbu- 
lent rough-wall  boundary  layer  suddenly  exposed  to  a smooth 

wall.  Computations  simulate  an  experiment  by  Antonia  and 
26 

Luxton,  for  which  an  incompressible  fluid  flows  past  a 
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Figure  16.  Comparison  of  computed  and  measured  effect  of 
surface  cooling  on  supersonic  boundary  layer 
transition;  Me  = 2.7,  Reoo=8 . 04  • 106/ft . 


flat  plate  perfectly  smooth  from  its  leading  edge,  x=0,  up 
to  x = 8 ft.  Between  x = 8 ft  and  x = 12  ft,  the  plate  is  covered 
with  distributed  (square)  roughness  elements  of  height  k - 1/8  in 
peak-to-peak  distance  between  the  elements  is  1/2  in.  Beyond 
x = 12  ft  the  plate  is  perfectly  smooth. 

The  objective  of  the  Antonia-Luxton  experiment  is  to  study 
boundary-layer  response  in  the  region  x > 12  ft.  The  boundary 
layer  undergoes  strong  departure  from  equilibrium  as  the  sur- 
face changes  from  rough  to  smooth.  In  particular,  Antonia  and 
Luxton  find  that  even  at  x=15.8  ft,  the  farthest  downstream 
location  at  which  they  measure,  the  boundary  layer  still  has 
not  achieved  equilibrium.  This  flow  hence  presents  a rigorous 
test  of  the  model. 

The  computation  was  performed  for  the  lowest  of  the  two-unit 
Reynolds  numbers  considered  by  Antonia  and  Luxton,  viz, 

U /v=l.ll*105  ft-1,  where  U is  freestream  velocity.  Cal- 
culation  was  initiated  at  x = 12  ft  from  the  measured  u,  <v,2>, 
and  <-u 'v' > profiles  . Figures  17-20  show  results  of  the  compu- 
tation (note  that  6S=2.8  in.  is  boundary-layer  thickness  at 
x = x = 12  ft).  Computed  velocity  profiles  virtually  duplicate 
measured  profiles  at  all  stations,  except  very  near  the  wall 
where  semilog  plots  indicate  10$  differences.  Computed  and 
measured  shape  factor  agree  to  within  3 $;  predicted  e^,  is 
within  20$  of  Preston  tube  measurements. 

The  most  significant  result  is  the  predicted  slow  approach  to 
equilibrium.  Inspection  of  computed  H and  c^.  show  that, 
consistent  with  measurements,  the  boundary  layer  is  still  out 
of  equilibrium  20  boundary  layer  thicknesses  downstream  of  the 
rough-wall/smooth-wall  juncture.  Discrepancies  between  com- 
puted and  measured  c^.  are  unalarming  because  Antonia  and  Luxton 

indicate  that  the  measurements  are  unreliable  for  (x-x„)/<5  <2. 
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Figure  17.  Comparison  of  computed  and  measured  velocity 
profiles  for  a turbulent  boundary  layer 
subjected  to  a sudden  change  in  surface 
roughness;  Uoo6s/v  = 2 . 6 • 10  4 . 
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Figure  19.  Comparison  of  computed  and  measured 

shape  factor  for  a turbulent  boundary 
layer  subjected  to  a sudden  change 
in  surface  roughness;  Uoo6s/v  = 2 . 6 • 10  4 . 
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Figure  20.  Comparison  of  computed  and  measured 

skin  friction  for  a turbulent  boundary 
layer  subjected  to  a sudden  change 
in  surface  roughness;  U 6 /v  = 2 . 6 • 10 4 . 
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However,  since  a rapid  initial  drop  in  cf  is  generally 
observed  for  such  a flow,  the  model's  predicting  this  subtle 
phenomenon  is  encouraging. 
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4.  BLUNT-BODY  APPLICATIONS 


The  preceding  section  presents  tests  of  the  turbulence/transi- 
tion model  for  well-documented  flows,  including  several  effects 
pertinent  to  transition  on  re-entry  vehicles.  Although  the 
tests  are  for  geometries  quite  different  from  re-entry  vehicles, 
the  tests  (because  the  flows  are  well-documented)  provide  a 
definitive  measure  of  model  accuracy.  However,  because  the 
primary  objective  of  this  study  is  to  devise  an  accurate  method 
for  predicting  transition  on  re-entry  vehicle  nosetips,  addi- 
tional tests  are  needed  which  include  the  special  conditions 
prevailing  for  flow  past  re-entry  vehicle  geometries.  In  this 
section,  we  hence  focus  upon  re-entry  vehicle  transition.  The 

section  consists  of  two  parts.  In  the  first  part  (Subsection  4.1), 
27 

the  PANT  ground-test  experiments  on  roughened  spheres  are 
simulated;  experimental  data  are  included  for  comparison.  The 
second  part  presents  results  of  a series  of  computations  test- 
ing the  model's  ability  to  predict  transition  sensitivity  to 
surface  roughness,  cooling  and  mass  addition;  high  Mach  number; 
freestream  noise  and  unit  Reynolds  number;  and  nonspherical 
geometry. 

4.1  PANT  GROUND- TEST  EXPERIMENTS 

In  our  first  re-entry  vehicle  application,  we  simulate  the  PANT 
experiments  in  which  a Mach  5 airstream  flows  past  a sphere- 
cone  of  nose  radius,  rN.  All  computations  are  for 

rN  = 2.5  inches  (59) 

while  the  ratio  of  surface  temperature,  T , to  freestream  total 
temperature,  T.  , is 

too 

Tw/Tt  = °*4  (60) 
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Boundary- layer-edge  pressure  is  given  by  the  modified  Newtonian 
pressure  distribution;  molecular  viscosity,  y,  is  determined 
from  the  Sutherland  law.  Freestream  turbulence  intensity  is 
assumed  to  range  from  about  0.10?  to  0.30?  so  that,  consonant 
with  the  arguments  presented  in  Appendix  A,  the  turbulent  mix- 
ing energy  at  the  boundary-layer  edge  is  given  by 


T* 


200 

3 


jf7 

Uavg 


2? 


(61) 


The  turbulence  scale  at  the  boundary-layer  edge  is 


£e  = .09  (62) 

Roughness  height,  k,  varies  from  0.1  mil  to  100  mils,  a much 

ii  5 

wider  range  than  considered  in  our  prior  applications.  ’ 


Figure  21  compares  computed  and  measured  transition  location 

P 8 

with  PANT  measurements;  the  PANT  correlation  is  also  shown 

for  reference.  Consistent  with  the  data,  model  predictions 

show  that  the  PANT. curve  provides  a reasonably  good  correlation 

for  (kT  /0T)  between  about  1 and  10.  Similar  to  the  measured 
e w t 

trend,  the  computed  Reg  asymptotes  to  a constant  value  for 

t 

both  small  and  large  roughness  heights.  Figure  22  shows  com- 
puted incipient  transition  unit  Reynolds  number,  Re^ , as  a 
function  of  k;  by  definition,  Re^  is  the  unit  Reynolds  number 
below  which  transition  occurs  downstream  of  the  sonic  point. 

Computed  values  of  Re  are  close  to  corresponding  measured 
27  1 

values . 


4.2  OTHER  TEST  CASES 

To  test  the  model's  ability  to  predict  transition  sensitivity 
to  the  various  phenomena  affecting  nosetip  transition,  further 
test  computations  have  been  made.  Experimental  data  sources 
were  unavailable  to  the  authors  at  the  time  the  computations 
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f 


the  PANT  wind  tunnel  experiments. 


k (mils) 


Figure  22.  Computed  incipient  transition  unit  Reynolds 
number:  M =5,  T/T.  =0.4. 


were  performed.  Table  2 summarizes  the  various  test  computa- 
tions; with  the  exception  of  Case  8 (see  Appendix  B)  , all  of 
the  test  cases  are  for  spherical  bodies.  For  each  case,  a 
transition  prediction  has  been  made;  then,  the  change  in  the 
value  of  the  "primary  variable"  required  to  change  the  transi- 
tion state  (i.e.,  laminar  to  turbulent  or  turbulent  to  laminar) 
has  been  determined  (provided  such  a change  is  possible). 

Subsections  4.2.1  through  4.2.6  present  results  of  the  compu- 
tations; the  concluding  subsection  summarizes  the  results. 

With  the  exception  of  freestream  turbulence  intensity  and  scale 
and  pressure  distributions , all  input  data  are  as  defined  in 
Table  2 and  Appendix  B.  All  cases  use  the  modified  Newtonian 
pressure  distribution  while  the  boundary-layer-edge  value  of 
the  turbulence  scale,  £.,  is  given  by  Equation  (62).  The  free- 
stream turbulence  intensity  varies  from  case  to  case.  Small 
values  (■v.Ol?)  are  used  for  flight-test  cases.  Larger  values 
( ''-I % ) as  defined  in  Figure  B1  of  Appendix  B are  used  for 
ground-test  cases. 

4.2.1  Combined  Effects  of  Roughness,  Cooling  and 
Mass  Addition 

We  begin  with  Test  Cases  1 through  3 which  are  clearly  ground- 
test  experiments.  The  working  fluid  is  air  with  constant 
specific  heat  ratio  y = 1.4;  viscosity  is  given  by  Sutherland's 
law.  The  only  assumed  input  datum  is  boundary-layer-edge 
turbulence  Intensity.  From  Figure  B1  of  Appendix  B,  the  value 
of  100/<u^2 >/Ue  is  found  to  be  0.75/5,  wherefore  we  assume 

T’  = 0.75 % (63) 

where  Uavp.  is  defined  in  Appendix  A. 

Table  3 summarizes  results  of  the  computations.  As  shown,  for 
the  baseline  case  in  which  k=5  mils,  B’  = 0.6  and  T /Th  =0.6, 

W ^ on 
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Table  2.  Blunt-Body  Test  Cases^" 


Case 

Number 

Effect  Analyzed 

Input  Conditions 

Primary 

Variable 

1 

Combined  effects  of 
surface  roughness, 
wall  cooling  and 
mass  addition 

Ptoo=  220  psia,  Ttoo=  850°R 

Moo  = ^ > rN  = ?" 

k=5  mils,  B'  = 0.6 

Tw/Ttoo  = 0 . 6 

k 

2 

Same  as  1 

Same  as  1 

B* 

3 

Same  as  1 

Same  as  1 

Tw/Ttoo 

4 

High  Mach  number 
combined  effects  of 
surface  roughness, 
wall  cooling  and 
mass  addition 

Flight  trajectory  in 
Table  B1  of  Appendix  B. 
Altitude  = 50  kft 
rN  = 0.75",  k = 0.4  mil 

Altitude 

5 

Same  as  4 

Same  as  4 except 
k = 15  mils,  non-sandgrain 
roughness  with  L/k  = 5 

Altitude 

6 

Freestream  noise  and 
unit  Reynolds  num- 
ber effects 

pt.=  260  psia,  Tfc  = 850°R 

r„-7" 

k = 5 mils  , B ' = 0 

Tw/Tto>  = 1 

Ptoo  and 
freestream 
turbulence 
intensity 
in  Fig- 
ure B1  of 
Appendix  B. 

7 

Surface  mass  addi- 
tion effects 

ptoo=  220  psia,  Ttoo  = 850°R 
M00=  6,  rN=  7” 
k = 0,  B'  * 0.6 
Tw/Ttoo  = 1 

Bf 

8 

Nonspherical  shape 
effects 

Same  as  1,  except  laminar 
stable  shape  as  given  in 
Table  B2  of  Appendix  B. 

k 

9 

Roughness  effects 
for  small  k/8 

Poo=  0.7  atm,  T « 495°R 
Mro=  3-14,  rN  = 6 . 5" 
k = .01  mil,  B f = 0 
Tw/Ttoo  = 0. 6 

k 

10 

Same  as  9 

Poo=  0.77  atm,  T = 517°R 
Moo=  2.73,  rN=  4ff 
k = .02  mil,  B'  = 0 
Tw^too  =0*75 

k 

In  all  but  Case  8,  the  body  is  a sphere;  B'  is  the  blowing  param- 
eter defined  by  B*  = pwvw/( pelIeCH)where  C^  is  Stanton  number;  in 
Case  5,  L is  peak-to-peak  distance  between  roughness  elements. 
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transition  occurs  very  close  to  the  stagnation  point,  viz,  at 
an  azimuthal  angle  of  13.8°.  The  predicted  transition 
coordinates  are 


Reft  = 87  (64) 

9t 

= 3-10  (65) 

et 

28 

where  is  Anderson's  modified  temperature  ratio  defined  by 

\p  = B'/10  + (1  + B'/4)T  /T  (66) 

W 6 

These  coordinates  are  very  close  to  the  values  predicted  by  the 
PANT  correlation  as  modified  by  Anderson  for  effects  of  mass 
addition. 


Table  3*  Summary  of  Computation  for 
Test  Cases  1 Through  3 


Test 

Case 

No. 

k(mils ) 

B 

! 

rp 

W 

Tt 

_ 

st 

rN 

♦t<°) 

M 

et 

Re 

9t 

T 

k , w 

0.  't 

t e 

k , . 

m 

5 

# 

6 

0 

6 

.240 

13.8 

0.30 

87 

3.87 

3-10 

■ 

■ 

l 

• 503 

28.8 

0.73 

193 

1.29 

1.04 

■ 

■ 

L 

.600 

34.4 

0.80 

207 

0.62 

0.50 

1 

■ 

1/2 

.6  00 

34.4 

0.80 

207 

0.31 

0.25 

l 

l 

0 

f 

r 

.612 

35.1 

0.82 

211 

0.00 

0.00 

2 

5 

0 

0 

6 

0.33 

84 

4.40 

4.40 

■ 

a 

I 

1 

im 

0 

8 

. 348 

19.9 

0.44 

115 

EEB 

2.55 

1 

1 

1 

* 

0 

9 

. 408 

23.4 

0.52 

128 

2.23 

■ 

■ 

■ 

1 

I 

■ 

1 

0 

.444 

25.4 

0.52 

153 

1.80 

■ 

a 

» 

1 

1 

1 

5 

.600 

34.4 

0.80 

132 

1.64 

1 

1 

1 

■ 

1 

1 

2 

0 

• 756 

43.3 

1.06 

99 

H 

1.68 

I 
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The  objective  of  the  Test  Cases  is  to  individually  vary  k,  B', 
and  T /T.  and  to  determine  if  variation  of  these  parameters 

W L go 

can  prevent  transition  from  occurring  ahead  of  the  sonic  point. 
Varying  k,  B'  and  T /T.  leads  to  the  following  conclusions: 

oo 

1.  It  is  not  possible  to  change  the  state  of  the  boundary 
layer  from  turbulent  to  laminar  by  decreasing  rough- 
ness height.  Even  in  the  perfectly  smooth  case  (k=0), 
transition  occurs  just  ahead  of  the  sonic  point. 

2.  Transition  occurs  for  all  blowing  rates,  including 
the  unblown  case  (B'=0). 

3.  T /T.  must  be  increased  to  2.0  to  cause  transition 

w r 

to  occur  beyond  the  sonic  point. 

4.2.2  High  Mach  Number  Combined  Effects  of  Roughness, 

Cooling  and  Mass  Addition 

Turning  now  to  Test  Cases  4 and  5,  we  test  the  model  under 
flight-test  conditions.  In  the  computations,  we  take  account 
of  high-temperature  properties  of  air  by  assuming  the  working 
fluid  to  be  a perfect  gas  with  constant  specific  heat  ratio 
y*1.2;  viscosity  is  given  by 

U = 2.5l8*icf10  T + 4.193*10'7  lb  *sec/ft2  (67) 

where  T is  temperature  in  degrees  Rankine.  All  other  input 
are  as  defined  in  Appendix  B. 

Two  freestream  turbulence  intensities,  T'  = .01?  and  .03?  are 

9 OO 

considered.  Table  4 summarizes  results  of  the  computations, 
including  inferred  values  of  the  boundary-layer-edge  intensity 
(Appendix  A) 

For  a roughness  height  of  0.4  mils  (Test  Case  4)  and  for 
turbulence  intensities  of  .01?  and  .03?,  no  transition  occurs 
at  50.2  kft.  Decreasing  altitude  fails  to  induce  transition 
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Table  Summary  of  Computations  for 
Test  Cases  ^ and  5 


even  at  37-7  kft  when  freestream  intensity  is  .01%.  However, 
when  .03%  freestream  turbulence  intensity  is  used,  transition 
first  occurs  at  ^0.7  kft. 


For  k = 15  mils  (Test  Case  5),  computations  have  been  performed 
two  ways.  First,  the  computations  make  no  special  allowance 
for  the  fact  that  the  roughness  is  non-sandgrain ; then,  compu- 
tations have  been  performed  with  surface-generated  turbulence 
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fluctuations  suppressed  as  a first  approximation  to  simulating 
distributed  roughness;  to  do  this,  we  replace  the  first  of 
Equations  (51)  by 

e = 0 at  y = 0 (68) 

The  following  results  have  been  obtained: 

No  Allowance  for  Non-Sandgraln  Roughness 

Transition  occurs  at  50.2  kft  for  freestream  intensities 
of  .01$  and  .03$.  For  both  intensities,  altitude  must 
be  increased  to  84.7  kft  in  order  to  change  the  state 
of  the  boundary  layer,  i.e.,  to  cause  transition  to 
occur  beyond  the  sonic  point. 

Surface  Disturbances  Suppressed 

When  freestream  intensity  is  .01$  the  flow  is  laminar 

at  50.2  kft;  for  .03$  freestream  turbulence,  transition 

occurs  at  50.2  kft.  Thus  for  T*  = .01$  we  decrease  alti- 

00 

tude  to  see  if  the  boundary  layer  will  transist;  however, 
the  flow  remains  laminar  all  the  way  down  to  37*7  kft. 

By  contrast,  we  increase  altitude  for  T^ = .03$  to  determine 
the  altitude  at  which  the  boundary  layer  remains  laminar 
ahead  of  the  sonic  point;  inspection  of  Table  4 shows 
that  no  transition  occurs  at  altitudes  greater  than 
60.2  kft. 

29 

Measurements  indicate  that  the  vehicles  considered  in  Test 
Cases  4 and  5 transist  at  nearly  the  same  altitude.  Although 
the  Test  Case  5 vehicle  transists  at  a higher  altitude  than 
that  of  Test  Case  4,  the  discrepancy  has  been  greatly  reduced 
by  using  Equation  (68).  Hence,  suppressing  surface-generated 
disturbances  appears  appropriate  for  analytically  simulating 
non-sandgrain  roughness. 
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4.2.3  Freestream  Noise  and  Unit  Reynolds  Number 

For  Test  Case  6,  computations  have  been  performed  for  the  five 
values  of  p.  at  which  /<u ' 2 >/U  is  given  in  Figure  B1  of 

^OG  6 0 

Appendix  B,  viz,  260,  210,  170,  130,  and  90  psia.  As  indi- 
cated in  Table  2,  this  test  case  has  neither  mass  addition 
nor  heat  transfer;  by  varying  pt  , we  vary  only  the  freestream 
unit  Reynolds  number.  Table  5 presents  results  of  the  compu- 
tations; as  shown,  transition  occurs  at  <J>t  = 21.3°  for  the 
baseline  case,  p = 260  psia;  total  pressure  must  be  reduced 
to  130  psia  to  achieve  laminar  flow. 


Table  5-  Summary  of  Computations  for  Test  Case  6 


Pt  (psia) 

<30 

10"6Re  (ft-1) 

OO 

T’  (%) 

st 

rN 

vo) 

M 

et 

Refl 

9t 

k /w 
9t  Te 

260 

5.47 

0.91 

• 372 

21.3 

0.47 

112 

3.37 

210 

4.06 

0.83 

. 516 

29.6 

0.67 

133 

2.59 

170 

3.29 

0.76 

. 720 

41.3 

0.99 

156 

1.94 

130 

2.52 

0.79 

◄ 

Laminar 

90 

1.74 

0.88 

◄ 



Laminar 

I 



► 

L_ 

4.2.4  Surface  Mass  Addition 

For  Test  Case  7>  we  isolate  the  effect  of  surface  mass  addi- 
tion on  nosetip  transition  under  ground-test  conditions.  As 
in  the  preceding  test  case,  we  have  no  heat  transfer;  in 
addition,  the  nosetip  is  perfectly  smooth  (k=0).  From 
Figure  B1  of  Appendix  B,  the  boundary-layer-edge  turbulence 
Intensity  is  0.75 %•  For  the  baseline  case  (B’  = 0.6),  transi- 
tion occurs  very  close  to  the  sonic  point  (see  Table  6). 
Decreasing  B'  causes  only  a slight  downstream  movement  of 


the  transition  point,  at 
(B'  = 0)  case. 


•mptoting  to  M =1.40  in  the  nonblown 


61 


Table  6.  Summary  of  Computations  for  Test  Case  7 




st 

rN 

*t(°) 

-p 

<D 

£ 

Re« 

9t 

msm 

.732 

41.9 

1.02 

232 

HI 

.780 

44.7 

1.10 

226 

Bfifl 

.900 

51.6 

1.35 

213 

■9 

.924 

52.9 

1.40 

205 

Bl 

.924 

52.9 

1.40 

205 

4.2.5  Nonspherical  Shape 

To  obtain  a measure  of  model-predicted  effects  of  body  geometry. 
Test  Case  8 considers  the  laminar-stable-shape  configuration 
(Appendix  B).  All  flow  conditions  match  those  of  the  Test 
Case  1 computations.  For  k=5  mils,  transition  occurs  at 
st/rN=  .280,  a bit  further  downstream  than  for  the  spherical 
body  considered  in  Test  Case  1.  Table  7 summarizes  results 
obtained  as  roughness  height  decreases  to  zero;  Figure  23  com- 
pares predicted  transition  location  for  the  laminar  stable  shape 
With  that  obtained  for  a sphere.  Decreasing  the  roughness 
height  to  k = 1 mil  is  sufficient  to  move  transition  to  a point 
Just  downstream  of  the  sonic  point;  further  decrease  of  k causes 
no  change  in  transition  point  location. 


Table  7-  Summary  of  Computations  for  Test  Case  8 


k(mils ) 

St 

rN 

M 

et 

Re 

6t 

k /w 
9t  Te 

-*-/* 

5 

.280 

U 

107 

3-65 

2.93 

1 

.693 

373 

0.35 

0.28 

1/2 

.693 

1.08 

352 

0.18 

0.15 

1/200 

.693 

1.08 

339 

0.0019 

0.0015 

0 

.693 

1.08 

339 

0 
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Figure  23-  Comparison  of  computed  effects  of  surface 
roughness  on  transition  location  for  the 
laminar  stable  shape  and  a sphere. 


4.2.6  Roughness  Effects  for  Small  k/0 

The  final  round  of  applications  (Test  Cases  9 and  10)  tests 
the  model's  accuracy  for  nosetip  transition  in  the  limit  of 
small  roughness  height.  Both  test  cases  are  for  vehicles  in 
free  flight;  Test  Case  9 simulates  the  experiments  of  Buglia 
while  Test  Case  10  corresponds  to  the  measurements  of  Hall, 
Speegle  and  Piland.^  Because  the  Mach  numbers  for  both  cases 

fall  in  the  supersonic  (rather  than  the  high  hypersonic)  range, 

the  working  fluid  is  assumed  to  be  air  with  y=  1.4  and  y given 
by  Sutherland's  law.  Freestream  turbulence  level  is  assumed 
to  be  .01%  so  that  according  to  Appendix  A, 

( . 088?  , Test  Case  9 

T'  = { .084%  , Test  Case  10  (69) 

For  Test  Case  9,  a roughness  height  of  .01  mil  causes  transi- 
tion Just  ahead  of  the  sonic  point  (see  Table  8).  At  transi- 
tion, the  predicted  s^/r^  is  .698;  by  comparison,  Buglia's 
measurements  indicate  that  sfc/r^  lies  between  0.66  and  0.78. 
Decreasing  k causes  only  a slight  movement  of  the  transition 
point;  for  the  perfectly  smooth  case,  transition  occurs  closer 
to,  but  still  upstream  of  the  sonic  point. 


Table  8.  Summary  of  Computations  for 
Test  Cases  9 and  10 


Test 

Case 

No. 

k(mils ) 

It 

rN 

V°) 

\ 

Refi 

9t 

. T 
k / w 

e/T 
t e 

.010 

£=- 

O 

o 

0.94 

654 

.0218 

1 

1 

.003 

40.7 

0.96 

664 

.0065 

■ 

■ 

.001 

• 711 

40.7 

0.96 

665 

.0022 

ft 

1 

0 

• 711 

40.7 

0.96 

666 

0 

.020 

- 

> 65 

— 

— 

— 

: 

.295 

1.128 

64.6 

1.839 

667 

.326 

■ 

1 

. 300 

1.044 

59.8 

1.627 

552 

.421 
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For  Test  Case  10,  the  flow  is  laminar  up  to  65°  from  the 
centerline  when  k=  .02  mil.  We  are  unable  to  integrate  beyond 
<(>  = 65°  because  the  Newtonian  pressure  distribution  fails  to 
provide  an  accurate  approximation  for  larger  values  of  $. 

Since  Hall,  Speegle  and  Piland  measure  transition  at  an  azimu- 
thal angle  between  about  75°  and  90°,  we  are  thus  unable  to 
make  a definitive  prediction  for  Test  Case  10.  However,  some 
insight  can  be  gained  by  using  larger  roughness  heights.  Com- 
putations for  k=  .295  mil  and  k * .300  mil  indicate  that  transi- 
tion may  occur  with  k=  .02  mil  at  a point  downstream  of  0 = 65°. 
That  is,  transition  occurs  with 


T 

■ w 


k . w 

e.  t 

t e 
, T 

= .326 

, ReQ  = 667 

for 

k = .295 

mil  ) 

(70) 

k / w 

e/T 
t e 

= .421 

, ReQ^  = 552 

for 

k = .300 

mil  ) 

These  points  are  close  to  the  PANT  correlation.  However,  we 
know  the  model  predicts  a finite  asymptotic  value  for  Re0  as 
k/8  -*-0.  This  asymptote  thus  occurs  for  Re.  >667,  and  hence 
for  <J>t>65°.  To  predict  the  transition  point,  the  computation 
would  have  to  be  repeated  using  a more  accurate  pressure 
distribution. 


4.2.7  Summary 

Based  on  the  results  obtained  for  the  PANT  simulations  of 
Subsection  4.1  and  for  the  ten  test  cases,  several  general 
observations  can  be  made. 


First , over  a very  wide  range  of  roughness  scale  heights, 
increasing  k destabilizes  a nosetip  boundary  layer.  However, 
the  change  in  the  boundary  layer's  stability  characteristics 
is  most  pronounced  for  a narrow  range  of  roughness  heights, 
viz,  0 . 5 <,  -Q- /\\>  < 10 . In  this  range,  Regt  varies  by  an  order 


of  magnitude.  Outside  this  range,  transition  is  much  less 
sensitive  to  variations  in  roughness.  For  example,  as  illus- 
trated in  the  PANT  simulations  and  in  Test  Cases  1,  9 and  10, 
when  roughness  height  decreases  to  very  small  values,  the 
transition  point  eventually  becomes  insensitive  to  roughness 
and  Re0  is  independent  of  k;  in  the  limit  k/0  -*■  0 transition 
location  depends  strongly  upon  freestream  turbulence  inten- 
sity. Similarly,  for  very  large  roughness  height,  Re0t  and 
st/rN  asymptote  to  constant  values  signifying  the  approach 
to  a completely  roughness  dominated  regime;  freestream  dis- 
turbances are  unimportant  in  this  regime.  Computations  for 
Test  Case  5 show  that  transition  induced  by  sandgrain  rough- 
ness is  quite  different  from  that  induced  by  non-sandgrain , 
distributed  type  roughness.  The  primary  difference  in  physi- 
cal mechanisms  active  for  these  two  kinds  of  roughness  is  in 
the  magnitude  of  surface  fluctuations  generated  by  the  rough- 
ness elements.  For  sandgrain  roughness,  local  flow  separations 
give  rise  to  strong  fluctuations  near  the  surface;  by  con- 
trast, for  distributed,  widely-spaced  roughness  elements,  sepa- 
ration will  be  more  orderly  and  less  likely  to  cause  random, 
turbulence-like  fluctuations  near  the  surface. 

Second,  the  strong  curvature  on  the  laminar  stable  shape  tends 
to  make  the  velocity  profile  more  stable.  As  shown  in  Fig- 
ure 23,  Reet  is  consistently  higher  at  all  roughness  heights 
for  the  laminar  stable  shape  than  for  a sphere. 

Third , as  illustrated  in  Test  Case  3,  cooling  destabilizes  a 
rough-wall  boundary  layer.  The  destabilization  occurs  because 
cooling  thins  the  boundary  layer  wherefore  k/0  increases. 

This  increase  in  k/0  makes  the  surface  appear  rougher  so  that 
Refl  decreases.  Test  Case  2 shows  that  surface  mass  addition 
has  a realtively  weak  effect  on  a rough-wall  boundary  layer 
under  ground-test  conditions.  The  computations  fail  to 
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Illustrate  an  interesting  possibility,  viz,  that  combined 
blowing  and  cooling  may  cancel  each  other's  effect  on  boundary- 
layer  stability.  That  is,  because  blowing  thickens  the  bound- 
ary layer  while  cooling  has  a thinning  effect,  a unique  com- 
bination of  blowing  and  cooling  may  permit  a constant  k/0  and 
consequently  a fixed  transition  point. 


Fourth , all  of  the  first  three  points  should  be  valid  for  the 
freestream  noise  environment  specified  in  Figure  B1  of  Appen- 
dix B.  In  fact,  this  relatively  noisy  environment  is  respon- 
sible for  the  relative  insensitivity  of  ReQ^_  to  roughness  and 
mass  addition  observed  in  Test  Cases  1 and  2.  In  both  cases 
the  boundary  layer  could  not  be  stabilized  by  any  change  in 

k or  B’.  The  limiting  k/0  0 value  of  Refl  , a limiting  value 

t 

strongly  controlled  by  freestream  noise,  is  achieved  at  finite 
values  of  k/0  and  B';  thus,  the  effectiveness  of  roughness  and 
mass  addition  is  Inhibited  by  freestream  noise. 


Fifth , the  first  four  conclusions  are  expected  to  hold  for  the 
range  of  freestream  Mach  numbers  from  2 to  20.  However, 
results  of  Test  Cases  4 and  5 suggest  that  nontrivial  depar- 
tures from  the  PANT  correlation  may  occur  at  Mach  20. 
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5.  DISCUSSION 


Results  presented  in  Sections  3 and  4 demonstrate  the  value 
of  the  turbulence-model  transition-prediction  method.  While 
the  program  objective  has  been  to  devise  an  accurate  and 
efficient  method  for  predicting  transition  on  re-entry  vehicle 
nosetips,  we  have  been  successful  in  predicting  transition 
for  other  aerodynamic  geometries  as  well.  The  model's  repre- 
senting more  than  one  class  of  flows  provides  hope  that  a 
universally  applicable  transition-prediction  method  has  been 
devised . 

Results  for  the  smooth-wall,  incompressible  FPBL  illustrate 
an  encouraging  aspect  of  the  model.  With  no  transition- 
specific  modifications  to  the  basic  turbulence  model  aside 
from  the  argument  used  to  set  the  value  of  the  closure  coeffi- 
cient X [Equation  (23)],  many  features  of  FPBL  transition  are 
accurately  predicted  including  transition  sensitivity  to  free- 
stream  turbulence  intensity,  transition  width,  skin-friction 
overshoot,  and  transitional  velocity  profiles.  In  addition, 
aside  from  boundary-condition  modifications  introduced  for 
numerical  reasons,  the  model  naturally  simulates  transition 
sensitivity  to  surface  roughness  and  surface  mass  addition. 

A primary  limitation  of  the  theory  is  its  inability  to  simulate 
spectral  effects.  The  closure  coefficients  in  the  model  equa- 
tions are  essentially  correlation  coefficients  which  have  been 
integrated  over  the  turbulent  spectrum.  Hence,  as  discussed 
in  Subsection  2.3  we  must  introduce  transition-specific  modi- 
fications in  order  to  accurately  predict  transition  sensitivity 
to  phenomena  which  reduce  the  range  of  unstable  frequencies 
such  as  pressure  gradient,  suction,  and  surface  heat  transfer. 
Nevertheless,  because  of  their  close  tie  to  linear-stability 
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theory,  the  modifications  proposed  in  Subsection  2.3  enhance 
the  model's  credibility  relative  to  earlier  versions  of  the 
model . 

Specific  areas  in  which  further  model  development  and  valida- 
tion would  be  beneficial  include  the  following:  (a)  surface 

disturbances  emanating  from  injected  mass,  (b)  differences 
between  sandgrain  and  non-sandgrain  roughness,  and  (c)  turbu- 
lence spectrum.  In  devising  the  surface  boundary  conditions  for 
surfaces  with  mass  addition,  we  have  assumed  the  Injected  gas 
to  be  free  of  turbulent  fluctuations.  This  is  probably  inaccu- 
rate for  materials  with  high  ablation  rates,  and  some  account 
of  surface  disturbances  should  be  made  in  such  applications. 
Similarly,  surface  disturbances  generated  by  roughness  ele- 
ments depend  upon  peak-to-peak  distance,  L,  between  roughness 
elements  as  well  as  upon  peak-to-valley  distance,  k.  Our 
rough-surface  boundary  conditions  are  designed  primarily  for 
sandgrain  roughness  which  is  the  special  limiting  case  in 
which  L/k^Od).  Because  many  materials  ablate  with  a range 
of  L to  k ratios,  further  development  of  the  rough-surface 
boundary  conditions  is  needed.  Finally,  the  model  masks  many 
of  the  subtle  wave  amplification  and  distortion  processes 
occurring  during  transition,  processes  which  are  strongly 
spectrum  dependent.  The  model's  credibility  would  be  enhanced 
considerably  if  some  direct  account  of  spectral  effects  were 
included . 

In  conclusion,  within  the  bounds  of  the  limitations  noted 
above,  the  model  is  ready  for  general  engineering  applications 
including  flight-test  cases.  Based  on  results  obtained  in 
this  two-ani-one-half  year  project,  we  expect  that  future 
applications  will  substantiate  this  claim. 
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APPENDIX  A 


I 


BOUNDARY-LAYER-EDGE  TURBULENCE  INTENSITY 
FOR  RE-ENTRY  VEHICLE  NOSETIPS 


To  determine  the  turbulence  intensity  at  the  boundary-layer 
edge,  T' , for  re-entry  vehicle  nosetips,  we  must  account  for 
two  effects: 

1.  Turbulence  amplification  near  the  stagnation 
streamline  caused  by  the  uniform  straining  region 
between  the  freestream  and  the  boundary-layer  edge. 

2.  Variations  in  T*  attending  rapid  acceleration  of 
the  fluid  away  from  the  stagnation  point. 

The  first  effect  is  pertinent  when  only  the  freestream  inten- 
sity, T^,  is  known.  The  second  effect  requires  special  atten- 
tion here  because,  for  simplicity,  we  use  a constant  average 
value  for  T'  in  all  computations;  accounting  for  the  second 
effect  permits  evaluation  of  an  appropriate  average  value. 

This  appendix  presents  the  method  used  to  account  for  these 
two  effects. 


In  all  of  the  blunt-body  computations  of  Section  4,  TJ^  is 
defined  in  terms  of  the  streamwise  fluctuating  velocity,  u^, 
as  follows: 

T'  = 100/<u'  2>'/U„  (Al) 

OO  OO  00 


Since  the  flow  is  normal  to  the  body,  the  definition  of  e 
[Equation  (12)]  implies  that 


T' 


200 

3 U 

OO 


( A2 ) 
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Similarly,  T'  and  eQ  are  related  as  follows 

e*5 

m,  - 200  ee 

3 Ue 


Therefore , 


(A3) 


(AU) 


32 

As  shown  by  Traci  and  Wilcox,  the  ratio  of  eg  to  em  on  the 

stagnation  streamline  primarily  depends  on  Reynolds  number. 

Figure  A1  shows  the  predicted  variation  of  ee/eoo  as  a function 

of  Reynolds  number  based  on  nose  radius,  Re_  . 

rN 


As  noted  above,  it  it,  convenient  to  assume  T’  is  constant 
throughout  the  nosetip  region.  This  is  not  rigorously  cor- 
rect, however,  since  Ug  varies  with  x while  eg  probably  changes 
only  slightly  as  x increases.  Thus,  some  average  value  of  T' 
must  be  devised  in  order  to  achieve  a meaningful  prediction. 

The  average  value  used  in  this  study  is  given  by 


T' 
T 1 


(A5) 


where  Uav^  is  the  average  velocity  between  the  stagnation  point 
(x=0)  and  the  sonic  point  (x=x*),  defined  by 

fX* 

U „ = ~w  U dx  A Iu«  ( A6) 

avg  x*  J e 2 e 

o 

where  U#  is  velocity  at  the  sonic  point.  The  integral  in 
Equation  (A^)  has  been  evaluated  by  noting  that  Ue  increases 
approximately  linearly  with  x.  Finally,  to  relate  U*  and  U , 
we  use  the  fact  that  total  temperature  is  constant  along  the 
t undary-layer  edge  wherefore 
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Figure  Al.  Stagnation  region  amplification  of  turbulent  mixing  energy  as 
function  of  nose-radius  Reynolds  number. 


APPENDIX  B 


INPUT  DATA  FOR  BLUNT-BODY  TEST  CASES 


Figure  B1  and  Tables  B1  and  B2  define  pertinent  input  data  for 
the  blunt-body  test  cases  of  Subsection  4.2.  The  quantity  B* 
is  the  blowing  parameter  defined  by 


B’ 


pwvw 

peUeCH 


where  is  Stanton  number. 


Shock  layer  disturbance  levels . 


I 

I 


Table  Bl.  Flight  Trajectory  Variables  for  Cases  4 and  5 


Altitude 

(ft) 

Velocity 
( ft/sec ) 

Stagnation 
Enthalpy 
( BTU/lb ) 

Freestream 

Density 

(lb/ft3 ) 
— 

pt 

C2 

(atm) 

B' 

T 

w 

(°R) 

V»t. 

98901 

22880 

10456 

.00105 

8.07 

li 

7082 

.45 

95338 

22851 

10430 

.00130 

9-97 

B il 

7142 

.46 

91778 

22814 

10397 

.00154 

11.8 

H !■ 

7191 

.47 

88222 

22771 

10357 

.00186 

14.2 

K 1 

7244 

.48 

84678 

22718 

10309 

.00221 

16.7 

.72 

7290 

.48 

81138 

22653 

10250 

.00260 

19-6 

.72 

7335 

.48 

77611 

22573 

10178 

.00310 

23.2 

.72 

7383 

.49 

74098 

22476 

10091 

.00368 

27.3 

.72 

7428 

.49 

70597 

22359 

9986 

.00448 

32.9 

.72 

7480 

• 50 

67119 

22213 

9856 

.00519 

37.6 

.72 

7516 

.50 

63668 

22037 

9700 

.00611 

43.6 

.71 

7555 

• 51 

60243 

21821 

9511 

.00719 

50.3 

.70 

7591 

.52 

56857 

21550 

9276 

.00846 

57-7 

.69 

7623 

.52 

53517 

21234 

9006 

.00986 

65-3 

.67 

7650 

.53 

50226 

20879 

8708 

.0116 

74.3 

.65 

7676 

.54 

46998 

20485 

8382 

.0135 

83.2 

.63 

7696 

.55 

43829 

20049 

8024 

.0157 

92.7 

.60 

7712 

• 56 

40732 

19572 

7652 

.0182 

102 

• 57 

7721 

.57 

37713 

19056 

7254 

.0210 

112 

.54 

7727 

• 58 

To  a good  approximation,  B'  may  be  considered  to  be  constant  over 
the  spherical  nosetip  at  a given  altitude. 
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Table  B2.  Laminar  Stable  Shape  Configuration 


Running 
Surface 
Distance , 

s* 

Local 

Body 

Angle, 

Local 

Body 

Curvature , 
d$/ds 

Local 
Radial 
Coordinate , 
r 

Local 

Axial 

Coordinate , 
z 

0.0 

90. 

-1. 

0.0 

0.0 

0.025217 

88.555 

-1.0002 

0.025214 

0.0003179 

0.050049 

87.132 

-1.0009 

0.050028 

0.0012524 

0.075137 

85.692 

-1.0022 

0.075066 

0.0028224 

0.10023 

84.250 

-1.0039 

0.10006 

0.0050215 

0.12506 

82.820 

-1.0061 

0.12473 

0.0078172 

0.15015 

81.372 

-1.0089 

0.14958 

0.011267 

0.17523 

79.920 

-1.0123 

0.17433 

0.0153^4 

0.20007 

78.477 

-1.0162 

0.19872 

0.019998 

0.22515 

77-013 

-1.0209 

0.22324 

0.025323 

0.25024 

75-541 

-1.0263 

0.24761 

0.031274 

0.27507 

74.077 

-1.0323 

0.27158 

0.037780 

0.30016 

72.588 

-1.0393 

0.29561 

0.044975 

0.32525 

71.089 

-1.0473 

0.31945 

0.052795 

0.35008 

69.592 

-1.0562 

0.34283 

0.061148 

0.37517 

68.067 

-1.0665 

0.36622 

0.070208 

0.40000 

66.541 

-1.0780 

0.38913 

0.079788 

0.42535 

64.967 

-1.0915 

0.41224 

0.090195 

0.45043 

63.387 

-1.1067 

0.43482 

0.10112 

0.47501 

61.817 

-1.1238 

0.45664 

0.11243 

0.50010 

60.188 

-1.1439 

0.47859 

0.12459 

0.52519 

58.527 

-1.1674 

0.50017 

0.13738 

0.55027 

56.830 

-1.1950 

0.52137 

0.15079 

0.57690 

54.980 

-1.2303 

0.54342 

0.16571 

0.60147 

53.221 

-1.2698 

0.56333 

0.18012 

0.62605 

51.400 

-1.3182 

0.58278 

0.19514 

0.65063 

49.502 

-1.3794 

0.60173 

0.21079 

0.67520 

47.505 

-1.4591 

0.62014 

0.22707 

0.70183 

45.191 

-1.5791 

0.63940 

0.24544 

0.72538 

42.959 

-1.7359 

0.65579 

0.26235 

0.75098 

40.220 

-2.0227 

0.67279 

0.28149 

0.77504 

37.085 

-2.6103 

0.68784 

0.30026 

0.80000 

31.955 

-4.4794 

0.70209 

0.32074 

0.82503 

27.348 

-2.2723 

0.71437 

0.34254 

0.85005 

24.777 

-1.4293 

0.72532 

0.36504 

0.87501 

23.047 

-1.0306 

0.73542 

0.38786 

0-90003 

21.746 

-0.80111 

0.74494 

0.41100 

0.92506 

20.710 

-0.65297 

0.75400 

0.43433 

0.95002 

19.854 

-0.54986 

0.76264 

0.45775 

0.97504 

19.123 

-0.47370 

0.77099 

0.48134 

1.0 

18.488 

-0.41545 

0.77903 

0.50497 

1.0250 

17.927 

-0.36928 

0.78685 

0.52874 

1.0501 

17.424 

-0.33179 

0.79446 

0.55264 

1.0751 

16.972 

-0.30098 

0.80184 

0.57649 

1.1000 

16.560 

-0.27511 

0.80904 

0.60039 

1.1251 

16.181 

-0.25297 

0.81611 

0:62446 

1.1501 

15.833 

-0.23401 

0.82299 

0.64845 

1.1750 

15.511 

-0.21751 

0.82973 

0.67248 

1.2000 

15.210 

-0.20303 

0.83634 

0.69655 

1.2251 

14.928 

-0.19016 

0.84287 

0.72078 

All  length  variables  have  been  normalized  by  the  stagnation  point 
radius  of  curvature  which  is  7.68  in.* 
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